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group data. We present results for the density of states, the density profile and 
the linear conductance. Two-particle backscattering leads to striking effects, which 
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the impurity and of the conductance at low energy scales is often preceded by a 
pronounced increase, and the asymptotic power laws are modified by logarithmic 
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I. INTRODUCTION 

One-dimensional metallic electron systems are always strongly affected by interactions. 
At low energy scales many observables obey anomalous power laws, known as Luttinger- 
liquid behavior, which is very different from conventional Fermi-liquid behavior describing 
most higher dimensional metals. 1,2 For spin-rotation invariant systems all power-law ex- 
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ponents can be expressed in terms of a single nonuniversal parameter K p . For Luttinger 
liquids with repulsive interactions [K p < 1) already one static impurity has a strong ef- 
fect at low energy scales, even when the impurity potential is relatively weak.— The 
asymptotic low-energy properties of Luttinger liquids with a single impurity have been 
investigated already in the 1990s.i^ For electron systems (spin-| fermions) with K p < 1 
the essential properties can be summarized as follows. The backscattering amplitude 
generated by a weak impurity is a relevant perturbation which grows as 

A (K p -l)/2 for a 

decreasing energy scale A. On the other hand, the tunneling amplitude through a weak 
link between two otherwise separate wires scales to zero as A Qs , with the boundary ex- 
ponent o>b = {Kp 1 — l)/2. At low energy scales any impurity thus effectively "cuts" 
the system in two parts with open boundary conditions at the end points, and physical 
observables are controlled by the open chain fixed point .» In particular, the local density 
of states near the impurity is suppressed as p ~ \u\ aB for \u\ — > 0, and the conductance 
vanishes as G(T) ~ T 2aB at low temperatures. We note that these power laws are strictly 
valid only in the absence of two-particle backscattering. In general they are modified by 
logarithmic corrections. The asymptotic behavior is universal in the sense that the expo- 
nents depend only on the properties of the bulk system, via K p , while they do not depend 
on the impurity strength or shape, except in special cases such as resonant scattering at 
double barriers, which require fine-tuning of parameters. 

The progress in the fabrication of artificial low- dimensional structures stimulated ad- 
vanced experimental verification of the theoretical predictions.— In an appropriate temper- 
ature and energy range Luttinger-liquid behavior can be expected in several systems with 
a predominantly one-dimensional character, such as organic conductors like the Bechgaard 
salts, artificial quantum wires in semiconductor heterostructures or on surface substrates, 
carbon nanotubes, and fractional quantum Hall fluids for chiral Luttinger liquids. For a 
correct interpretation of experimental data it would be helpful to have theoretical input 
beyond asymptotic power laws, which are valid only at sufficiently low energy scales. It 
is not always clear whether the asymptotic Luttinger-liquid behavior is well developed 
before finite size effects and interactions with the three-dimensional environment become 
important. 
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Recently, a functional renormalization group (fRG) method has been developed for 
a direct treatment of microscopic models of interacting fermions with impurities in one 
dimensio n n i 12 i 13 which not only captures correctly the universal low-energy behavior, but 
allows one to compute observables on all energy scales, yielding thus also nonuniversal 
properties, and in particular an answer to the important question below which scale the 
asymptotic power laws are actually valid. The method has been applied to the spinless 
fermion model with nearest-neighbor interaction on a one-dimensional lattice, supple- 
mented by various types of impurity potentials. The most relevant observables such as 
the local density of states , 11 ! 12 the density profile,— and the linear conductanc e) 13 ! 14 were 
calculated. The truncation of the fRG hierarchy of flow equations employed in these 
works is valid only for sufficiently weak interactions. However, a comparison with exact 
numerical results from the density matrix renormalization group (DMRG) 15 showed that 
the truncated flow equations are generally rather accurate also for sizable interaction pa- 
rameters. The fRG captures complex crossover phenomena at intermediate scales, such 
as the temperature dependence of the conductance through a resonant double barrier . 13 i 16 
It can also be applied to other (than chain) geometries, such as mesoscopic rings threaded 
by a magnetic flu?*^ or Y junctions.— 

In this work we extend the fRG method for interacting Fermi systems with a single or 
few impurities to spin-| fermions and apply it to the extended one-dimensional Hubbard 
model. For fermions with spin, vertex renormalization is crucial to take into account that 
two-particle backscattering of fermions with opposite spins at opposite Fermi points scales 
to zero in the low-energy limit. By contrast, for spinless fermions the effects of a single 
static impurity are captured qualitatively already within the lowest order truncation of 
the fRG hierarchy of flow equations, where the renormalized vertex is approximated by 
the bare interaction. Two-particle backscattering scales to zero only logarithmically, and 
thus gives rise to logarithmic corrections to the asymptotic power laws. 

The analysis of fixed-point models, which yields the ultimate low-energy behavior, 
predicts a power-law decay of the local density of states near an impurity or boundary of 
systems with K p < 1. However, a non-selfconsistent Hartree-Fock and a DMRG study of 
the density of states near the end of a finite Hubbard model chain with open boundary 
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conditions revealed that the power-law suppression at the lowest scales can be preceded by 
a pronounced increase of spectral weight . 19 ' 20 A similar crossover behavior can therefore 
also be expected for the density of states near an impurity, at least for a sufficiently 
strong one, as we indeed obtain in this work from the fRG. For the conductance, a 
renormalization group analysis of the g-ology model by Matveev et al.— showed that two- 
particle backscattering can lead to an increase as a function of decreasing temperature 
before the asymptotic suppression sets in. 

The paper is organized as follows. In Sec. II we introduce the microscopic model and 
derive the corresponding fRG flow equations. In Sec. Ill we present results for spectral 
properties of single-particle excitations near an impurity or boundary, the density profile, 
and transport properties. We conclude with a summary and an outlook in Sec. IV. 

II. MODEL AND FLOW EQUATIONS 

A. Microscopic model 

As a microscopic model for the bulk electron system we choose the one-dimensional 
extended Hubbard model with a nearest-neighbor hopping amplitude t, a local interaction 
U, and a nearest-neighbor interaction U' . The bulk system is supplemented by a site or 
hopping impurity. The total Hamiltonian is given by 

H = -* X) ( C 5+V c i* + C ]<x c i+V ) + U Yl n ^ n ^ + U ' X n i n i +1 + Him ? ' W 
j> j j 

where cj CT and c J(7 are creation and annihilation operators for fermions with spin projection 

a on site j, while rij a = c^c^, and rij = + rij\ is the local density operator. For 

the (nonextended) Hubbard model the nearest-neighbor interaction vanishes. A local site 

impurity on site jo is modeled by H imp = V nj , and a hopping impurity by the nonlocal 

potential H imp = (t - 1') J2 a ( c j +i,a c j ,a + c ) ,a S'o+V ) > sucn tnat tne no PP in g amplitude 
t is replaced by t' on the bond linking the sites jo and jo + 1 • In the following we will set 
the bulk hopping amplitude t equal to one, that is all energies are expressed in units of t. 

In the absence of impurities, the Hubbard model can be solved exactly using the Bethe- 
ansatz,— while the extended Hubbard model is not integrable. The Hubbard model is a 
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Luttinger liquid for arbitrary repulsive interactions at all particle densities except half- 
filling, where the system becomes a Mott insulator: 1,2 The phase diagram of the extended 
Hubbard model is more complex. Away from half-filling, it is a Luttinger liquid at least 
for sufficiently weak repulsive interactions.— For the Hubbard model the Luttinger-liquid 
parameter K p can be computed exactly from the Bethe ansatz solution.— 

For the calculation of transport properties a finite interacting chain (with sites 1, . . . , L) 
is coupled to noninteracting leads at both ends. The influence of the leads on the inter- 
acting chain can be taken into account by incorporating a dynamical boundary potential 



VJ- (*«.) = + (l - ^ - J -±-^ J {SlJ + Slj) , (2) 

in the bare propagator Gq of the interacting chain.— The parameter /i is the chemical 
potential, which is related to the density n in the leads by fio = —2 cos kp with kp = mr/2. 
Uncontrolled conductance drops due to scattering at the contacts between leads and the 
interacting part of the chain can be avoided by switching off the interaction potential 
smoothly near the contacts. In addition, interaction induced bulk shifts of the density 
have to be compensated by a suitable bulk potential.— 



B. Flow equations 

We now extend the fRG scheme derived and used previously for spinless Fermi systems 
with impuritie o llll2lls to electrons, that is spin-| fermions. We make use of equations and 



12 and Ref. 
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without repeating 



procedures described already in detail in the articles Ref. 
the derivations here. 

We use the one-particle irreducible (1PI) version of the fRG . 23 ' 24 ^ The starting point 
is an exact hierarchy of differential flow equations for the 1PI vertex functions, which is 
obtained by introducing an infrared cutoff A in the free propagator and differentiating the 
effective action with respect to A. Since translation invariance is spoiled by the impurity, 
we use a Matsubara frequency cutoff instead of a cutoff on momenta. The cutoff is 
sharp at T = 0— and smooth for T > . 13 The hierarchy is truncated by neglecting 
the contribution of the three-particle vertex to the flow of the two-particle vertex. The 
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coupled system of flow equations for the two-particle vertex T and the self-energy S is 
then closed. The contribution of the three-particle vertex to T A is small as long as T A is 
sufficiently small. 

We neglect the influence of the impurity on the flow of the two-particle vertex, such 
that T A remains translation invariant. While this is sufficient for capturing the effects 
of isolated impurities in otherwise pure systems, it is known that impurity contributions 
to vertex renormalization become important in macroscopically disordered systems.— We 
also neglect the feedback of the bulk self-energy into the flow of P\ which yields only a 
very small correction at weak coupling. The two-particle vertex is parametrized approxi- 
mately by a renormalized static short-range interaction^ in order to reduce the number of 
variables in the flow, which would be unmanageably large otherwise. This approximation 
is exact at the beginning of the flow and fully captures the nonirrelevant parts of the 
vertex in the low-energy limit. The self-energy generated by the simplified vertex is then 
static (frequency independent) and its spatial dependence can be treated fully, that is 
without resorting to another simplified parametrization. Transport properties are com- 
puted by coupling the interacting model to noninteracting leads as described in Ref. [kI 
The conductance is obtained directly from the one-particle Green function, since current 
vertex corrections vanish in our approximation for T A . 

We now describe the parametrization of the spatial (or momentum) dependences of the 
two-particle vertex T A for spin-| fermions, employing a natural extension of our previous 
parametrization for the spinless case.— We consider only spin-rotation invariant lattice 
systems with local and nearest-neighbor interactions. This includes the extended Hubbard 
model. 

For a spin-rotation invariant system the spin structure of the two-particle vertex can 
be decomposed into a singlet and a triplet part: 



r 



.A 




(3) 



with 
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°i ,°- 2 ;cti,o- 2 



(4) 
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Since the total vertex is antisymmetric in the incoming and outgoing particles, the singlet 
part has to be symmetric and the triplet part antisymmetric. 

Proceeding in analogy to the case of spinless fermions,— we first list momentum com- 
ponents of the vertex with all momenta at ±kp. For the triplet vertex the antisymmetry 
allows only one such component 

9t = ^t\k F ,-k F ;k F ,~k F ■ (5) 

For the singlet vertex there are several distinct components at ±fcp. Since we will ne- 
glect the influence of the impurity on the vertex renormalization, the renormalized vertex 
remains translation invariant. Hence the momentum components are restricted by mo- 
mentum conservation: k[ + k' 2 = k\ + &2, modulo integer multiples of 2tc. The remaining 
independent (not related by obvious symmetries) components are 

o A = r A 

js2 s\k F ,—k F ;k F ,—k F ' 



9s4 F ' s \kp,kp;kp,kp i (6) 

and in the case of half-filling, for which kp = 7r/2, also 

9s3 = r s | 7r /2, 7 r/2;-7r/2,-7r/2 • (?) 

The labels 2,3,4 are chosen in analogy to the conventional g-ology notation for one- 
dimensional Fermi systems.— In order to parametrize the vertex in a uniform way in all 
cases, we will include the umklapp component g^ 3 not only at half-filling, but at any 
density. The effect on the other components is negligible for the range of interactions and 
fillings considered. 

Extending our treatment of the spinless case,— we now parametrize the vertex by 
renormalized local and nearest-neighbor interactions in real space. For the triplet part, 
there is no local component, and only one nearest-neighbor component compatible with 
the antisymmetry, namely 

U t = ^t\j,j+l;j,j+l ) ( 8 ) 

which has the same form as the nearest-neighbor interaction in the spinless case. Note 
that ^tljj+i-jj+i does not depend on j, and is equal to ~^t\jj-i-jj-v F° r the symmetric 
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singlet part, there is one local component 



s\j,j;j,j ' 

and three different components involving nearest neighbors 



U s T s | 7,7:7,7 > (9) 



u' s A 


-- r A 




= r A 


w A 


- r A 



(!°) 

For the Hubbard model, the bare vertex is purely local and the initial condition for the 
vertex is given by [7 A ° = 2U, while all the other components vanish. For the extended 
Hubbard model, U' S A ° = U' t A ° = U' are nonzero. 

The triplet vertex is parametrized by only one renormalized real space coupling, which 
leads to a momentum representation of the form 

?U^-mm = 2f/ t A t cos (^ - - cos (^ - fei)] 5<tS k2 ,K + k> , (11) 

where the Kronecker delta implements momentum conservation (modulo 27r). The flow- 
ing coupling U' t k is thus linked in a one-to-one correspondence to the Fermi momentum 
coupling (yf A by 

( 7 A = 2f/; A [l-cos(2A; F )] ) (12) 

as in the spinless case.— In the singlet channel we have found four real space couplings, 
that is one more than necessary to match the three singlet couplings in momentum space, 
9s2i 9s3i 9s4- We choose to discard the interaction W S A , because it does not appear in 
the bare Hubbard model, where it is generated only at third order in U, while the pair 
hopping P A appears already in second order perturbation theory. Fourier transforming 
the remaining interactions yields the singlet vertex in k-space 



s| k' 1 ,k' 2 ;ki,k2 



£/ A + 2U' S A [cos(X - fci) + cos(A;2 - h)] + 2P A cos(A; 1 + k 2 



U k 1 +k 2 ,k' 1 +k' 2 

(13) 
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from which we obtain a linear relation between the momentum space couplings g A 2 , g A 3 , 



g^ 4 and the renormalized interaction parameters [/ A , U' S A , P A : 

g A = U A + 2U' S A [1 + cos(2M] + 2P S A 
gf 3 = U A - AU' A - 2P A 

gi = U A + iU' A + 2P A cos(2A; F ) . (14) 

The determinant of this linear system is positive for all hp, except for kp = and tt. 
Hence the equations can be inverted for all densities except the trivial cases of an empty 
or completely filled band. 

We can now set up the flow equations for the four independent couplings U' t A , f/ A , 
U' S A , and P A which parametrize the vertex. Consider the case T = first. Inserting the 
spin structure (J3J) into the general flow equation for the two-particle vertex, Eq. (18) in 



Ref. [12J, and using the momentum representation for a translation invariant vertex, the 
flow equation for the singlet and triplet vertices T A , a — s,t, can be written as 



^rW^-sE £ /f( pp + ™ + PH ' 

(.1=4- A h h' — e t J 



(15) 



=±A b,b'=s,t ' 

with the particle-particle and particle-hole contributions 



— C^ h , Gp(iu) Gp +kl _ k ^ (iuj) r 6 | fe / jP+fcl _ fe /. fcl)p r 6 ,| Pife /. J , +fel _ A ./ 



k 2 



PH' — C™, Gp(iuj) G p+kl __ k ,^{iu]) ^ A \k' 2 ,p+k 1 -k' 2 M,P^b'\p,k' 1 -,p+ki-k' 2 ,k2 ■ (^) 

The coefficients C a ,w are obtained from the spin sums as 




3/4 (17) 
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Note that we have neglected the self-energy feedback in the flow of T , such that only 
bare propagators Go enter. On the right hand side of the flow equation we insert the 
parametrization (jllj) for Tf and ()13|) for T^. The flow of the triplet vertex T^ k , k ,. ki fca 
is evaluated only for (k[, k' 2 , k\, k 2 ) = (hp, —kp, kp, —kp) as in ©, which yields the flow 
of g£, while the flow of the singlet vertex T^ k , k ,. ki k2 is computed for the three choices 
of k[, k' 2 ,ki, k 2 which yield the flow of gy 2 , gf 3 , g^ A . Using the linear equations (fT2"j) and 
()14|) to replace the couplings g A by the renormalized real space interactions on the left 
hand side of the flow equations, we obtain a complete set of flow equations for the four 
renormalized interactions U[ , U A , U' S A , and P A of the form 

d ^ =EE h *A«>) U$ U A , (18) 

lu=±A a', a" 

where a = 1,2,3,4 labels the four different interactions. The functions h®, a ,j{uj) can be 
computed analytically by carrying out the momentum integrals in (I15|) via the residue 
theorem; the flow equations can then be solved numerically very easily. The expressions 
become too lengthy to be reported here; for details we refer the interested reader to 
the thesis by one of the authors. 27 For finite systems the momentum integral should be 
replaced by a discrete momentum sum; however, this leads only to negligible corrections 
for the physical observables presented in Sec. III. 

After computing the flow of the real space interactions, one can also calculate the flow 
of the momentum space couplings g A by using the linear relation between the two. In 
the low-energy limit (small A) one recovers the one-loop flow of the g-ology model, the 
general effective low-energy model for one-dimensional fermions.— In addition, our vertex 
renormalization captures also all nonuniversal second order contributions to the vertex at 
±kp from higher energy scales. 

In Fig. ^ we show results for the renormalized real space interactions together with the 
corresponding momentum space couplings, as obtained by integrating the flow equations 
for the Hubbard model at quarter-filling and T = 0. Note that the couplings converge to 
finite fixed-point values in the limit A — > 0, but the convergence is very slow, except for the 
momentum space couplings g A 3 and g A A . This can be traced back to the familiar behavior 
of the so-called backscattering coupling g K L = \ (g A 2 — g A ), that is the amplitude for the 
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exchange of two particles with opposite spin at opposite Fermi points. Backscattering is 
known to vanish logarithmically in the low-energy limit for spin-rotation invariant spin- 
| Luttinger liquids.- We emphasize that this logarithmic behavior is not promoted to 
a power law by higher order terms beyond our approximation. By contrast, the linear 
combination of couplings which determines the Luttinger-liquid parameter K p converges 
very quickly to a finite fixed-point value (see below). 

Due to the above parametrization of the vertex by real space interactions which do 
not extend beyond nearest neighbors on the lattice, the self-energy generated by the flow 
equations is frequency independent and tridiagonal in real space. Inserting the spin and 
real space structure of T A into the general flow equation for the self-energy, Eq. (16) in 



Ref. 



12l one obtains 

d „ A 1_ 
r 

w=±A r=±l 



^ A . = -- £ [u$Gf^)HK L + ut) E 
M = -h E [ W - + p " ] . (is) 

uj=±A 

where G A = (Gq 1 — E A ) -1 . These equations can be solved very efficiently," 8 so that very 
large systems with up to 10 7 sites can be treated. 

At T > the Matsubara frequencies are discrete, and a sharp frequency cutoff there- 
fore leads to discontinuities in the flow. To avoid ambiguities and numerical problems 
associated with these discontinuities one may choose a smooth frequency cutoff as in 



Ref. 



13l . Alternatively, one can rewrite the Matsubara sum as a frequency integral with 



a suitable weight function, as described in detail in the Appendix. The latter procedure 
leads to a particularly simple and numerically convenient extension of the flow equations 
to finite temperatures. In the flow equations (jIBJl and (JTTJJ) one has to replace uj = ±A 
by L) = ±c<; A , where uj a is the Matsubara frequency which is closest to A. The function 
Kx'a"{ u ) remains the same. At T > the flow equations can be solved for systems with 
up to 10 5 sites without extensive numerical effort. 

When using a sharp cutoff with G A defined above it makes no difference whether 
the bare impurity potential is put into Gq 1 or S A ; we choose to include it in the initial 
condition of E A at A = oo. Due to the slow decay of G at large frequencies, the integration 
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of the flow equation for X from A = oo to A = A yields a contribution which remains 
finite even in the limit A — > oo.- 12 For the extended Hubbard model this contribution is 
given by = U/2 + 2U' for j = 2, . . . , L-l and = X£° L = U/2 + U' . The numerical 
integration of the flow is started at a sufficiently large A with S A ° as initial condition. 



C. Calculation of K p 

The Luttinger-liquid parameter K p can be computed from the fixed-point couplings 
as obtained from the fRG. A relation between the fixed-point couplings and K p can be 
established via the exact solution of the fixed-point Hamiltonian of Luttinger liquids, the 
Luttinger model. Since the above simplified flow equations yield not only the correct 
low-energy asymptotics to second order in the renormalized interaction, but contain also 
all nonuniversal second order corrections at at higher energy scales, the resulting K p 
is obtained correctly to second order in the interaction. 

The Luttinger-liquid parameter K p is given by 



"" l/ l + (^4 + ^ 2 )/(^) ' { } 

The coupling constants g P 2 and g P 4 parametrize forward scattering interactions in the 
charge channel (which is spin symmetrized) between opposite and equal Fermi points, 
respectively. They are related to the bare singlet and triplet vertices of the Luttinger 
model by 

9p2 ^ ( Ts| kp,—kp;kp,—kp ~\~ ^ >r )t\ kp,—kp;kp ,—kp^) 

1 

9f>4 — -^ls\kp,kp;kF,kp ■ (21) 

These bare vertices are identical to the dynamical forward scattering limits of the full 
vertex T. On the other hand, the vertex T A obtained from the fRG with a frequency 
cutoff yields the static forward scattering limit for A — > 0. 12 For the Luttinger model, the 
static forward scattering limit of the vertex can be computed from the effective interac- 
tions D p2 {q,iv) and D p ^(q,iv), which are defined as the sum over all particle-hole chains 
with the bare interactions g P 2 and g P 4.— The summation becomes a simple geometric 
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series if one introduces symmetric and antisymmetric combinations g p ± = g p ^ ± g p i and 
D p ±(q, iv) = D pi (q, ii/)±D p2 (q, if). The static limit of the effective interaction D p ±(q, iv) 
yields the relation 

9 * P± = 1 - W(vr^) (22) 
between the Luttinger-model couplings g p ± and the fixed-point couplings 

s;± = ^k*4±U* 2 + 3^)] (23) 
from the fRG with frequency cutoff. Inverting (J22j) one obtains 



K ' = \U- a ;j(^)- (24) 

The Fermi velocity Vf can be computed from the self-energy for the translation invariant 
pure system as in the spinless case, 12 using the momentum representation of the flow 
equations (JTJJ). 

The results for K p from the above procedure are correct to second order in the bare 
interaction for the Hubbard model and also for the extended Hubbard model. While the 
flowing couplings g^ 2 and g£ converge only logarithmically to their fixed-point values for 
A — > 0, the linear combination g^ 2 + which enters K p converges much faster. 

In Fig. 121 we show results for K p for the Hubbard model as obtained from the fRG 
and, for comparison, from the exact Bethe ansatz solution.— The truncated fRG yields 
accurate results at weak coupling except for low densities and close to half-filling. In 
the latter case this failure is expected since umklapp scattering interactions renormalize 
toward strong coupling, even if the bare coupling is weak. At low densities already the bare 
dimensionless coupling U/vp is large for fixed finite U, simply because Vp is proportional 
to n for small n, such that neglected higher order terms become important. Note, for 
comparison, that for spinless fermions with a fixed nearest-neighbor interaction the bare 
dimensionless coupling at the Fermi level vanishes in the low-density limit. 

For the extended Hubbard model Fig. |3] shows a comparison of fRG results for K p to 
DMRG data. 29 The fRG results are exact to second order in the interaction and are thus 
very accurate for weak U and U'. Results from a standard one- loop g-ology calculation 26 
deviate quite strongly already for U' > 0.5. In the g-ology approach interaction processes 
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are classified into backward scattering (gi±), forward scattering involving electrons from 
opposite Fermi points (g2±), from the same Fermi points (<?4_i_), and umklapp scattering 
(<?3_i_)- All further momentum dependences of the vertex are discarded. This is justified 
by the irrelevance of these momentum dependences in the low-energy limit, but leads to 
deviations from the exact flow at finite scales, and therefore to less accurate results for 
the fixed-point couplings. 

The flow of gf]_, i = 1, ...,4, is plotted in Fig. in the upper panel for the quarter- 
filled Hubbard model with bare interaction U — 1, and for the extended Hubbard model 
with U' = U I y/2 in the lower. The fRG result is compared to the result from a one-loop 
g-ology calculation. The backscattering coupling g\± vanishes logarithmically in both 
cases, as expected for the Luttinger-liquid fixed point.— For the pure Hubbard model 
the good agreement with g-ology results stems from the purely local interaction in real 
space, since in that case pronounced momentum dependences of the vertex develop only 
in the low-energy regime where the g-ology parametrization is a good approximation. 
By contrast, for the extended Hubbard model momentum dependences of the vertex 
which are not captured by the g-ology classification (except for small A) are obviously 
more important. A generalization of the g-ology parametrization of the vertex to higher 
dimensions, which amounts to neglecting the momentum dependence normal to the Fermi 
surface, is frequently used in one-loop fRG calculations in two dimensions.— The above 
comparison indicates that this parametrization works well for the pure Hubbard model, 
but could be improved for models with nonlocal interactions. The parametrization of the 
vertex by an effective short-range interaction used here could be easily extended to higher 
dimensions, where it will probably yield more accurate results, too. The relevance of an 
improved parametrization of the vertex beyond the conventional g-ology classification has 
also been demonstrated in a recent fRG analysis of the phase-diagram of the half-filled 
extended Hubbard model.— 
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III. RESULTS FOR OBSERVABLES 

In this section we present and discuss explicit results for the spectral properties of 
single-particle excitations, the density profile, and the conductance for the Hubbard and 
extended Hubbard model with a single impurity, as obtained from the solution of the 
fRG flow equations. We also analyze excitations and density oscillations near a boundary, 
which corresponds to an infinite site impurity or a vanishing weak link. A comparison 
with DMRG results is made for the spectral weight at the Fermi level and for the density 
profile. For details on the computation of the relevant observables from the solution of 



the flow equations we refer to Ref. 
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A. Single-particle excitations 

Integrating the flow equation for the self-energy E A down to A = yields the physical 
self-energy E and the single-particle propagator G = (Gq 1 — From the Green 

function G the properties of single-particle excitations can be extracted. We focus on the 
local spectral function given by 

pAuo) = -- ImGjjicu + i0 + ) . (25) 

7T 

For a finite system this function is a finite sum of 5-peaks of weight w\j, where A labels 
the eigenvalues of the effective single-particle Hamiltonian defined by G. Dividing the 
spectral weight w\j by the level spacing yields the local density of states Dj(u). Even- 
odd effects due to finite-size details in the spectral weight are averaged out by averaging 
over neighboring eigenvalues. 

For uo — > the spectral weights and the local density of states near a boundary or 
impurity are ultimately suppressed according to a power law with the boundary exponent 

a B = + 1 (26) 

2K P 2K a K } 

with K a = 1 for spin-rotation invariant systems.— However, due to the slow logarith- 
mic decrease of the two-particle backscattering amplitude, the fixed-point value of K a is 
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reached only logarithmically from above. Hence, we can expect that the asymptotic value 
of a.B is usually reached only very slowly from below. 

The local density of states at the boundary of a quarter-filled Hubbard chain, com- 
puted by the fRG, is shown in Fig.EJfor various values of the local interaction U. Contrary 
to the expected asymptotic power-law suppression, the spectral weight near the chemi- 
cal potential is strongly enhanced. The predicted suppression occurs only at very small 
energies for sufficiently large systems. In the main panel of Fig. El the crossover to the 
asymptotic behavior cannot be observed, as the finite size cutoff ~ rrvp/L is too large. 
Results for a larger system with L = 10 6 sites at U = 2 in the inset show the crossover 
to the asymptotic suppression, albeit only at very small energies. The dependence of the 
boundary spectral weight at the Fermi level on the system size L is plotted in Fig. |BJ 
The L-dependence of the spectral weight at zero energy is expected to display the same 
asymptotic power-law behavior for large L as the ^-dependence discussed above. Instead 
of decreasing with increasing L, the spectral weight increases even for rather large systems 
for small and moderate values of U. For U > 2 the crossover to a suppression is visible 
in Fig. |U1 For U = 0.5 only an increase is obtained up to the largest systems studied. 
The crossover scale depends sensitively on the interaction strength U; for small U it is 
exponentially large in Vp/U. 

The above behavior of the spectral weight and density of states near a boundary of 
the Hubbard chain, that is a pronounced increase preceding the asymptotic power-law 
suppression, is captured qualitatively even by the Hartree-Fock approximation . 19 i 20 This 
is at first sight surprising, as the Hartree-Fock theory does not capture any Luttinger- 
liquid features in the bulk of a translation invariant system. The initial increase of Dj(uS) 
near a boundary is actually obtained already within perturbation theory at first order in 
the interaction,— 



space interaction, and z the number of spin components. For spinless fermions (z — 1) 
with repulsive interactions the coefficient in front of the logarithm is always positive such 




(27) 



where D°-{uj) is the noninteracting density of states, V(q) the Fourier transform of the real 
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that the first order term leads to a suppression of Dj(u). For the Hubbard model, one 
has z = 2 and V"(0) — 2V(2kp) = —U is negative for repulsive U. Hence, at least for weak 
U the density of states increases for decreasing uo until terms beyond first order become 
important. For the extended Hubbard model, V'(O) - 2V{2k F ) = 2U'[1 - 2 cos(2A; F )] - U, 
which can be positive or negative for U, U' > 0, depending on the density and the relative 
strength of the two interaction parameters. At quarter-filling V(0) — 2V(2kp) is negative 
and therefore leads to an enhanced density of states for U' < U/2. 

Using g-ology notation, one can write V"(0) — 2V(2kp) = g2± — 2gij_, which reveals that 
substantial two-particle backscattering (gi± > g2±/2) is necessary to obtain an enhance- 
ment of Dj{uj) for repulsive interactions. Backscattering vanishes at the Luttinger-liquid 
fixed point, but only very slowly. In case of a negative V(0) — 2V(2kp) the crossover 
to a suppression of Dj(u) is due to higher order terms, which are expected to become 
important when the first order correction is of order one, that is for energies below the 
scale 

/ 2ttvf \ 

cu c = e^exp — ~ , (28) 

\V{0)-2V{2k F )J V ' 

corresponding to a system size L c = ttvf/uj c . The scale u c is exponentially small for weak 
interactions. A more accurate analytical estimate of the crossover scale from enhancement 
to suppression has been derived for the Hubbard model within Hartree-Fock approxima- 
tion in Ref. |2JJ. In a renormalization group treatment uj c is somewhat enhanced by the 
downward renormalization of backscattering. 

A comparison of fRG results with DMRG data 32 for the spectral weight at the Fermi 
level is shown in Fig.[7l for a boundary site in the upper panel, and near a hopping impurity 
of strength t' = 0.5 in the lower. The agreement improves at weaker coupling, as expected, 
and is generally better for the impurity case, compared to the boundary case. The larger 
errors in the boundary case are probably due to our approximate translation invariant 
parametrization of the two-particle vertex. Boundaries and to a minor extent impurities 
spoil the translation invariance of the two-particle vertex. Although the deviations from 
translation invariance of the vertex become irrelevant in the low-energy or long-distance 
limit, and therefore do not affect the asymptotic behavior, they are nevertheless present 
at intermediate scales. This feedback of impurities into the vertex increases of course 
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with the impurity strength and is thus particularly important near a boundary. The scale 
for the crossover from enhancement to suppression of spectral weight discussed above 
depends sensitively on effective interactions at intermediate scales and can therefore be 
shifted considerably even by relatively small errors in that regime. 

With the additional nearest-neighbor interaction in the extended Hubbard model it is 
possible to tune parameters such that the two-particle backscattering amplitude becomes 
negligible. In that case the asymptotic power-law suppression of spectral weight should be 
free from logarithmic corrections and accessible already for smaller systems and at higher 
energy scales. The bare backscattering interaction in the extended Hubbard model is 
given by gi± — U + 2U' cos(2kp) and therefore vanishes for U' = — U/[2 cos(2fcp)], which 
is repulsive for U > if n > 1/2. In a one-loop calculation a slightly different value of U' 
has to be chosen to obtain a negligible renormalized g^ ± for small finite A, since the flow 
generates backscattering terms at intermediate scales even if the bare g\± vanishes. In 
Fig.|H]we show fRG and DMRG results^ for the spectral weight of the extended Hubbard 
model at the Fermi level near a hopping impurity. In the upper panel a generic case with 
sizable backscattering is shown, while the parameters leading to the curves in the lower 
panel have been chosen such that the two-particle backscattering amplitude is negligible 
at low energy. Only in the latter case a pronounced suppression of spectral weight is 
reached already for intermediate system size, similar to the behavior obtained previously 
for spinless fermions with nearest-neighbor interaction . 11 ' 12 This is also reflected in the 
energy dependence of the local density of states near the impurity. For parameters leading 
to negligible two-particle backscattering as in Fig.|U]the suppression of the density of states 
sets in already at relatively high energies and is not preceded by any interaction-induced 
increase. Note also that the fRG results are much more accurate for small backscattering, 
as can be seen by comparing the agreement with DMRG data in the upper and lower 
panel of Fig. |S] especially for larger U. This indicates that the influence of the impurity on 
the vertex flow, which we have neglected, is more important in the presence of a sizable 
backscattering interaction. 

In the case of a negligible backscattering amplitude, the spectral weight at the Fermi 
level approaches a power law without logarithmic corrections for accessible system sizes 
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if the impurity is sufficiently strong. The power law is seen most clearly by plotting the 
effective exponent ot(L), that is the negative logarithmic derivative of the spectral weight 
with respect to the system size. Fig. [TU1 shows ot{L) on the site next to a site impurity of 
strength V for the extended Hubbard model with U — 1, U' — 0.65, and n = 3/4. The 
backscattering amplitude is very small for these parameters. The fRG results approach the 
expected universal ^-independent power law for large L, but only very slowly for small 
V. For a weak bare impurity potential V, the crossover to a strong effective impurity 
occurs only on a large length scale of order V 2 ^ Kp ~ 1 ^ .— For V = 0.1 this scale is obviously 
well above the largest system size reached in Fig. The Hartree-Fock approximation 
also yields power laws for large L, but the exponents depend on the impurity parameters. 
This failure of Hartree-Fock theory was already observed earlier for spinless fermions.— 

The effective exponent obtained from the fRG calculation agrees with the exact bound- 
ary exponent to linear order in the bare interaction, but not to quadratic order. To 
improve this, the frequency dependence of the two-particle vertex, which generates a fre- 
quency dependence of the self-energy, has to be taken into account. This is also necessary 
to describe inelastic processes and to capture the anomalous dimension of the bulk sys- 
tem. These effects could be included in an improved scheme by inserting the second 
order vertex into the flow equation for the self-energy without neglecting its frequency 
dependence. 

B. Density profile 

Boundaries and impurities induce a density profile with long-range Friedel oscillations, 
which are expected to decay as a power law with exponent (K p + K a )/2 at long distances, 
where K a = 1 for spin-rotation invariant systems.— For weak impurities linear response 
theory predicts a decay as \j — jo\ 1 ~ Kp ~ K ' 7 at intermediate distances. 

The density profile has to be computed from an additional flow equation, since the 
local density is a composite operator whose renormalization is not well described by the 
propagator obtained from the truncated flow of S A . The flow equation for can be 
derived by computing the shift of the grand canonical potential fl A generated by a small 
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field <pj coupled to the local density. Its general structure at T = is described in Ref. 
From Eqs. (38) and (40) in that article one can easily obtain the concrete flow equation 
for the case of the extended Hubbard model. 

As an additional benchmark for the fRG technique, we compare in Fig. ^2 fRG and 
DMRG results for the density profile rij for a quarter-filled Hubbard chain with L = 128 
lattice sites and open boundaries. Friedel oscillations emerge from both boundaries and 
interfere in the center of the chain. The fRG results have been shifted by a small constant 
amount to allow for a better comparison of the oscillations. Note that the mean value 
of rij in the tails of the oscillations deviates from the average density by a finite size 
correction of order 1/L, which is related to the asymmetry of the oscillations near the 
boundaries. 

The long-distance behavior of the density oscillations as obtained within the fRG 
scheme has been analyzed in detail for spinless fermions in Ref.tl2|. For fermions with spin, 
asymptotic power laws can be identified only for special parameters leading to negligible 
two-particle backscattering. In general, the asymptotic behavior of Friedel oscillations 
is realized only at very long distances, and the power laws are modified by logarithmic 
corrections. 



For the computation of the conductance a finite interacting chain is connected to two 
semi-infinite noninter acting leads, with a smooth decay of the interaction at the contacts. 
The presence of leads modifies the propagator in the interacting region only via the 
boundary potential V lead , Eq. (J2J). In linear response the conductance is given b y 13 ' 34 



with \t(e,T)\ 2 = [4 - (/i + e) 2 ] \G hL (e,T)\ 2 , and / the Fermi function. The factor two 
is due to the spin degeneracy. Within our approximation scheme, S A has no imaginary 
part, which implies that there are no vertex corrections, such that the conductance is fully 
determined by boundary matrix element of the single-particle propagator Gi^l{£, T).— 



C. 



Conductance 




(29) 
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For a system of spinless fermions with a single impurity it was already shown that the 
conductance obtained from the truncated fRG obeys the expected power laws, in partic- 
ular G(T) oc T 2olb at low T, and one-parameter scaling behavior . 13 i 14 The corresponding 
scaling function agrees remarkably well with an exact result for K p = 1/2, although the 
interaction required to obtain such a small K p is quite strong. The more complex temper- 
ature dependence of the conductance in the case of a double barrier at or near a resonance 
is also fully captured by the fRG . 13 i 16 

Fig. shows typical fRG results for the temperature dependence of the conductance 
for the extended Hubbard model with a single strong site impurity (V = 10). Similar 
results were obtained for a hopping impurity. The considered size L = 10 4 corresponds 
to interacting wires in the micrometer range, which is the typical size of quantum wires 
available for transport experiments. For U' — the conductance increases as a function 
of decreasing T down to the lowest temperatures in the plot. For increasing nearest- 
neighbor interactions U' a suppression of G(T) at low T becomes visible, but in all the 
data obtained at quarter-filling the suppression is much less pronounced than what one 
expects from the asymptotic power law with exponent 2a b- By contrast, the suppression 
is much stronger and follows the expected power law more closely if parameters are chosen 
such that two-particle backscattering becomes negligible at low T, as can be seen from 
the conductance curve for n = 3/4 and U' = 0.65 in Fig. The value of K p for these 
parameters almost coincides with the one for another parameter set in the plot, n = 1/2 
and U' = 0.75, but the behavior of G(T) is completely different. Note that at T ~ nvp/L 
finite size effects set in, as can be seen at the low T end of some of the curves in the figure. 
An enhancement of the conductance due to backscattering has been found already earlier 
in a renormalization group study of impurity scattering in the g-ology models 

Results for the conductance of the extended Hubbard model with a hopping impurity 
with various amplitudes t' are shown in Fig. The bulk parameters have been chosen 
such that the two-particle backscattering is practically zero at low T. From the plot of the 
logarithmic derivative of G(T) in the upper panel one can see that for a strong impurity 
(small t') the conductance follows a well defined power law G(T) oc T 2oiB over a large 
temperature range. For intermediate t' the curves approach the asymptotic exponent at 
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low T from below, but do not reach it before finite size effects lead to a saturation of 
G(T) for T < ttvf/L. For the weakest impurity in the plot, t' = 0.95, the conductance 
remains very close to the unitarity limit. However, the plot of the logarithmic derivative of 
l — G/(2e 2 /h) in the lower panel of Fig. fH?l shows that l — G/(2e 2 /h) increases as T Kp ~ l for 
decreasing T, as expected for a weak impurity in the perturbative regime.— The effective 
exponents indicated by the two horizontal lines in the figure deviate from the exact values 
(determined from the DMRG result^ for K p ) by about 20 % in the case of 2a# and only 
by 5 % for K p — 1. Results for the conductance of a wire with a double-barrier impurity 
at or near resonance will be presented in a forthcoming publication.— 

IV. CONCLUSION 

We have derived a fRG-based computation scheme for the one- dimensional extended 
Hubbard model with a single static impurity, extending previous work for spinless 
fermion o 12113 to spin-| fermions. The underlying approximations are devised for weak 
short-range interactions and arbitrary impurity potentials. Various observables have been 
computed: the local density of states near boundaries and impurities, the density profile, 
and the temperature dependence of the linear conductance. Results have been checked 
against DMRG data, for those observables and system sizes for which such data could 
be obtained. The general agreement is good at weak coupling, but for intermediate 
interaction strengths with sizable two-particle backscattering and strong impurities the 
deviations are significantly larger than for spinless fermions. We suspect that the ne- 
glected influence of impurities on vertex renormalization at high energy and short length 
scales is more important for fermions with spin. 

Two-particle backscattering of particles with opposite spin at opposite Fermi points 
leads to two important effects, not present in the case of spinless fermions. First, the 
expected decrease of spectral weight and of the conductance at low energy scales is of- 
ten preceded by an increase, which can be particularly pronounced for the density of 
states near an impurity or boundary as a function of uo. For the density of states near 
a boundary this effect has been found already earlier within a Hartree-Fock and DMRG 
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study of the Hubbard model , 19 > 20 and for the conductance by a renormalization group 
analysis of the g-ology model.— Second, the asymptotic low-energy power laws are usually 
modified by logarithmic corrections. In the extended Hubbard model the backscattering 
can be eliminated for a special fine-tuned choice of parameters. Then the results are very 
similar to those for spinless fermions. For weak and intermediate impurity strengths the 
asymptotic low-energy behavior is approached only at rather low scales, which are acces- 
sible only for very large systems. This slow convergence was observed already for spinless 
fermion a 11 ! 12 ! 13 and holds also in the absence of two-particle backscattering. 

For systems with long-range interactions backscattering is strongly reduced compared 
to forward scattering. This seems to be the case in carbon nanotubes.— Hence, the conduc- 
tance can be expected to follow the asymptotic power law at accessible temperature scales 
for sufficiently strong impurities in these systems, as is indicated also by experiments.— 
However, the effects due to two-particle backscattering should be observable in systems 
with a screened Coulomb interaction. 
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APPENDIX A: FREQUENCY CUTOFF AT FINITE TEMPERATURE 

In this appendix we derive a convenient implementation of a sharp frequency cutoff 
at finite temperature. The idea is roughly to rewrite the Matsubara sum as an integral 
over a piecewise constant function and then to introduce a sharp cutoff on this continuous 
frequency as usual. 

In the 1PI scheme^ the flow equation with frequency cutoff at zero temperature has 
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the form (assuming no frequency shift along the loop)^ 

a£- = J ^4(kl-A)/(e e (|^|-A),^,T A ). (Ai) 

Here T A is the generating functional for the 1PI vertex functions, and f(t,u, T A ) repre- 
sents the right-hand side of the flow equations including the momentum integral, but with 
the integral over the Matsubara frequency uj written explicitly. Q e (x) is a step function 
smoothed on a scale e, and S e (x) = d x Q e (x). The u integral gives finite contributions near 
|a>| = A, hence T A is a continuous function of A. Then the limit e — > can be performed 
using Morris' lemma,— 



~~8A 



/ %Ls(\u,\-A) [ 1 dtf(t,u;,r A ) = ±- fdtf(t^). 

J l ^ Jo Z7r ,»io 



(A2) 

u=±\ Jo 

The right-hand side is independent of the explicit cutoff function Q e (x), therefore the 
vertex functions in T A are smooth not only in A but also as functions of all external 
frequencies. 

In our previous work at finite temperature^ the flow equation (jAl|) with discrete 
Matsubara frequencies reads 
<9T A 



OA 



Tj2^(\oo n \ -A)/(6 e (K|-A),u; n ,r A ). (A3) 



In the limit e — > the right-hand side contains a 5 function and T A jumps as A passes 
u n , hence the Morris lemma which requires continuity of T A cannot be applied. We have 
therefore used a smooth cutoff, which renders the numerics relatively slow. Instead, we 
now propose a new sharp cutoff which allows to apply the Morris lemma. This reduces 
the runtime significantly and enables us to access another order of magnitude in system 
size and temperature. 

First, we rewrite the Matsubara sum as an integral over a continuous frequency uj with 
a weight function peaked in non-overlapping neighborhoods of width r\ around each u n , 



T ^2fM= / dujT^Sr, 

J n 



[u-u n )f(M(u)) (A4) 



where J dio 5 v (u — u n ) = 1 and M(u) returns the discrete Matsubara frequency u) n closest 
to uj. Hence, / is a piecewise constant function of a continuous variable uj. At this stage, 
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the cutoff function in u is introduced in the bare action. This leads to the flow equation 
— - = / duuTj2^-^n) S e (\u\ - A) f(e e (\u\ - A),MH,T A ). (A5) 

n 

If we take the limit 77 — ^ first we again obtain equation (|A3|) . which is the smooth cutoff 
limit. For finite 77, however, Morris' lemma gives 

J duTj2^-^n)6(\u\-A) ^ dt f(t, M(uj), T A ) 

= TVM Wn |-A) [ dtf(t,u n ,T A ). (A6) 

For each tu n separately this is an autonomous differential equation, hence the result is 
independent of the shape of S v (x). A convenient choice is a box of height l/(27rT) and 
width 2nT centered around the Matsubara frequency u n , that is, 5 v (x) = l/(2irT) for 
\x\ < 71T and otherwise, which leads to the final result 

<9T A 



i- r * t a ). 



9A 



(A7) 



Comparison with equation ()A2|1 shows that the only change necessary at finite temper- 
ature is to replace the loop frequency uo = ±A on the right-hand side by the nearest 
discrete Matsubara frequency. We have checked numerically that indeed this new sharp 
cutoff gives the same results as the previous smooth cutoff for the conductance G{T) 
curves in a dramatically reduced runtime. 
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FIG. 1: (Color online) Vertex flow for the Hubbard model at quarter-filling (n = 1/2) and 
U = 1; upper panel: flow of the renormalized real space interactions, lower panel: flow of the 
momentum space couplings. 
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FIG. 2: (Color online) Luttinger-liquid parameter K p for the Hubbard model as a function of 
electron density. Results from the fRG are compared to exact results from the Bethe ansatz. 
The upper curves are for U = 1 and the lower ones for U = 2. 
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FIG. 3: (Color online) Luttinger-liquid parameter K p for the quarter-filled extended Hubbard 
model as a function of U' for U = 0.5 and {7=1. Results from the fRG are compared to DMRG 
data and to results from a one-loop g-ology calculation. 
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FIG. 4: (Color online) Flow of the vertex on the Fermi points (in g-ology notation) at quarter- 
filling and U = 1; upper panel: Hubbard model, lower panel: extended Hubbard model with 
U' = U /y/2. The fRG flow is compared to the one-loop g-ology flow; note that in the upper 
panel fRG and g-ology results almost coincide. 
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FIG. 5: (Color online) Local density of states at the boundary of a Hubbard chain of length 
L = 4096 at quarter-filling and various interaction strengths U; the inset shows results for U = 2 
and L = 10 6 at very low to. 
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FIG. 6: (Color online) Spectral weight at the Fermi level at the boundary of a quarter-filled 
Hubbard chain as a function of system size L, for various different interaction strengths. 





32 



- 1 1 — i i i i 1 1 1 1 1 — i i i i 1 1 1 1 1 — i i i i 1 1 1 1 1 — i i i i 1 1 



e — e — o 




U = 1 



_i i I i i I i i I i i 



0.7 



0.6 



3 

^ 0.5 



0.4 - 



- 1 1 — I I I I I 1 1 1 1 — I I I I I 1 1 1 1 — I I I I I 1 1 1 1 — I I I I I I 




J I I I I I I I 



10 1 



10 2 



10 3 
L 



10 4 



10° 



FIG. 7: (Color online) Spectral weight at the Fermi level near a boundary {upper panel) and a 
hopping impurity t' = 0.5 (lower panel) as a function of system size L for the Hubbard model 
at quarter-filling and different interaction strengths U; results from the fRG (open symbols) are 
compared to DMRG data (filled symbols). 




FIG. 8: (Color online) Spectral weight at the Fermi level near a hopping impurity t' = 0.5 as a 
function of system size L for the extended Hubbard model with U' = U / v2, for various choices 
of U ; upper panel: n = 1/2 (leading to sizable backscattering) , lower panel: n = 3/4 (leading to 
small backscattering); results from the fRG (open symbols) are compared to DMRG data (filled 
symbols). 
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FIG. 9: (Color online) Local density of states near a hopping impurity t' = 0.5 in an ex- 
tended Hubbard model with density n = 3/4 and interaction U' = U/y/2 (leading to a small 
backscattering interaction) for various choices of U. The size of the chain is L = 4096. 




FIG. 10: (Color online) Logarithmic derivative of the spectral weight at the Fermi level on the 
site next to a site impurity of strength V in the center of the chain as a function of system size 
L, for the extended Hubbard model with interaction parameters U = 1, U' = 0.65, and density 
n = 3/4; here the filled symbols are fRG, the open symbols Hartree-Fock results. 



35 




FIG. 11: (Color online) Density profile rij for the Hubbard model with 128 sites and interaction 
strength U = 1 at quarter-filling; fRG results are compared to DMRG data. 




FIG. 12: (Color online) Temperature dependence of the conductance for the extended Hubbard 
model with L = 10 4 sites and a single site impurity of strength V = 10, for a Hubbard interaction 
U = 1 and various choices of U'\ the density is n = 1/2, except for the lowest curve, which has 
been obtained for n = 3/4 and U' = 0.65 (leading to a very small backscattering interaction); 
the dashed line is a power-law fit for the latter parameter set. 
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FIG. 13: (Color online) Logarithmic temperature derivative of the conductance (upper panel) 
and of its deviation from the unitarity limit (lower panel) for the extended Hubbard model with 
L = 10 4 sites and various hopping impurities. The density is n = 3/4, interaction parameters 
are U = 1 and U' = 0.65. The dashed horizontal lines highlight power-law behavior. 



